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CROSS-REFERENCES TO RELATED APPLICATIONS 
[01] This application claims the benefit of U.S. Provisional Application No. 
60/220,604, filed July 25, 2000, which disclosure is incorporated by reference. 



BACKGROUND OF THE INVENTION 
[02] Pitch measurement is the measurement of the distance between two 
similar features. This definition applies to both the physical specimen and its image, that is, 
the corresponding data set collected from the specimen using a metrology or imaging system. 
The sample pitch is the spacing of similar elements in a specimen such as a grating. In most 
metrology and imaging systems, when a specimen with known sample pitch is measured or 
imaged, the spacing between similar features is rendered in the data set or image and the 
pitch in such data set or image can be used to determine the scale in the data set, for example, 
the magnification of the image. Thus, pitch measurement can be used for metrology or 
imaging system calibration to set the scale. If the scale is known in advance of data 
collection, a measurement of the sample pitch in a specimen with otherwise unknown pitch 
can be obtained from the pitch in the data set collected from the sample. 

[03] Prior art methods for pitch measurement in metrology and imaging 
systems fall into two categories. In the first category of prior art, measurement methods 
determine the pitch in the data set in a manner similar to linewidth measurement. Hence, a 
discussion of linewidth measurement methods follows. 

[04] In most metrology and imaging systems the signal is formed as a result 
of interaction of excitations with a specimen or sample and the detection of all or part of the 
signal generated as a result of the interaction. The excitations, in some cases referred to as 
probes, can be mechanical probes, electrons, photons, ions, phonons or other forms of 
radiation. In magnetic resonance imaging, for example, the excitations are magnetic field 
pulses. The detected signal can consist of charged particles, photons, phonons, or other 
observables such as temperature or field amplitude. What is common among most metrology 
and imaging systems is that they all create one, two, or three dimensional spatial images of 
specimen. That is, the systems produce a vector function f of the spatial coordinates r in 
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some region of the space. The function f describing the signal can be a function of one, two, 
or three dimensional coordinates as f(x), f(x,y), or f(x,y,z). In addition, the function f can be 
a vector itself having several components. For example, a color image is represented as three 
functions over space: red, green and blue. These systems often produce information that is 
stored in computer memory in some digital format as a single column or multidimensional 
matrices or graphics formats for images, though the method of pitch measurement also 
applies to continuous signals such as photographs since they can be scanned and converted 
into digital forms. The data collected from metrology or imaging instruments will be referred 
to as a data set. This term includes gray scale images obtained from any number of systems 
such as CCD cameras and scanning electron microscopes. The one-dimensional data is 
referred to as a profile data. This term is intended to include scan lines collected with 
scanning probe microscopes or other scanning instruments capable of producing such data. 

[051 The features in a data set are rendered from features on the sample, 
though sometimes such data set features can result from noise in the scan or other electronics 
or from environmental factors such as vibration and other interferences. When data set 
features correspond to those of the sample, information about the physical features is 
obtained from their rendition in the collected data. In the ideal case, the feature as rendered 
in the data set would contain information that is identical to that of the physical feature. 
However, metrology and imaging systems have imperfections consisting of, but not limited 
to, sharpness of excitation (probe sharpness), finite interaction volume of excitation and 
sample, detection limits and inefficiencies, finite depth of focus, signal to noise ratio 
limitations, diffraction, and image distortion. In the data set, the portions of the data set 
formed from the sharp edges of physical sample features are smeared by these imperfections. 

[06] As an example, consider a specimen such as a diffraction grating, 
which consists of a series of parallel lines of equal width and equal spaces between the lines. 
Further assume that the metrology or imaging acquires data from a portion of the specimen 
consisting of parts of two or more lines. The linewidth measurement is the determination of 
the distance between the two edges of a physical feature, for example, the left edge and the 
right edge of a line in a diffraction grating. In the data set, linewidth or feature width 
measurement is essentially the distance between the corresponding two edges in the data set. 
Therefore, methods of measuring the width of a physical feature are based on methods that 
determine the position of two edges in the data set. With such metrology and imaging 
systems, the problem is in determining where the actual edge location lies along a complex 
waveform of the detected signal. Sometimes, the edge position in the data set is modeled 
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based on the physics of the system and the particular experimental conditions including 
probe-sample interactions. But more often, the edge position in the data set is assigned by 
applying any number of arbitrary edge detection algorithms to the data set or a representation 
of the data set. Common edge detection algorithms are threshold (absolute or percentage), 
5 maximum derivative or slope, second derivative, s curve fit, linear regression, peak-to-peak 
distance, and centroid-to-centroid distance. All these algorithms consider a transition in the 
detected signal (low-to-high or high-to-low transition, or the combination of the two) and 
assign to some point in the signal transition interval a pixel position which is designated as 
the edge location in the data set, and hence corresponds to the feature edge location in the 
10 sample. Generally these methods implement subpixel interpolation schemes. The net result 
of smearing and arbitrary edge assignment is that the linewidth measurement tends to be 
biased and suffers from lack of accuracy, 
yp [07] The methods in the first category of prior art for pitch measurement are 

n nearly identical to linewidth measurement as described above. The pitch in the data set is 

: 

jfi measured as the difference between the locations of two similar edges of the same type in the 

£n data set. That is, a left-edge to left-edge distance or a right-edge to right-edge distance, or 

i i \ 

^ ' peak-to-peak distance and so on. This type of prior art is depicted in Figure 1, where the 
« original data set (not shown) is reduced to a one-dimensional profile data 2. In this example, 
ffj transitions in the profile data represent the edges in the sample. Two similar edge regions 12 
30 and 14 (both left edges in this example) in the profile data are determined. The pixel 

positions 16 and 18 of the two similar edges 12 and 14 are considered. The positions 16 and 
18 are determined by application of any one of edge detection algorithms identical to the 
algorithms used for edge position determination in linewidth measurement listed before. The 
pitch in the data set is then calculated as the difference between the position of two similar 
25 edges in the data set, P = r2 - rl, where P is the measured pitch in the data set, rl is the pixel 
location 16 of the first edge as determined by a particular algorithm, and r2 is the pixel 
location 18 of the second edge as determined by the same algorithm. The edge pixel 
positions are often determined to a subpixel range of values by interpolation. In summary, 
the first category of prior art requires determination of two similar edge positions in the data 
30 set. 

[08] The smearing of the features in the image or data set affects linewidth 
and pitch measurements differently. In a linewidth measurement, the results depend on the 
magnitude of the broadening of the edge in the data set as caused by the particular physical 
effects and experimental conditions that cause the broadening of the edges at the time of data 
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collection. This fact is well known, and for this reason, image magnification calibration to 
the feature with a known width is not recommended and in modem metrology tools, 
instrument image magnification calibration is accomplished with measurement of pitch of a 
specimen with known physical pitch. Due to the fact that the smearing alters the two similar 
5 edges in the data set in a similar manner, the effects of broadening of the two edges tend to 
cancel out in pitch measurement. For this reason pitch measurements are referred to as self- 
compensating or unbiased. More specifically, such pitch measurements do not require a 
physical model for edge determination in order to be accurate and are far less sensitive to 
image imperfections arising from finite probe size (excitation-sample interaction volume). 
10 However, it is known in the prior art that all pitch measurement algorithms are not identical 
in performance. At issue is the undesired contribution of the measurement algorithm to the 
uncertainty of the pitch measurement. This contribution limits the precision or repeatability 

S with which the pitch can be measured, and consequently the precision with which the scale of 

M= the metrology or imaging system can be calibrated. 

U 5 [09] The National Institute of Science and Technology (NIST) has 

y 1 developed a version of the pitch measurement algorithm belonging to the first category of 
5 prior art for use in scanning electron microscopes. This algorithm is described in an article 
~1 entitled "A New Algorithm for the Measurement of Pitch in Metrology Instruments", N. F. 
rU Zhang et al, Proceedings of SPIE, Vol. 2725, pp. 147-158, 1996. In the NIST algorithm, the 
q20 two edge locations are determined using regression to fit two straight lines to the signal in the 
^ transition regions for each edge, but the method places an additional restriction on the slopes 
of the lines fitted to the two similar edges in the data. The fitted lines must have identical 
slopes in the regression. 

[10] The second category of prior art consists of methods that apply Fourier 
25 transform to the data set as disclosed in U.S. Patent No. 4,818,873 issued to Glen A. Harriot. 
These methods determine the location of the peak in the amplitude of the Fourier transform 
of the data set or detected signal. Since the peak in Fourier Transform of the signal 
corresponds to the fundamental spatial frequency in the image or data set, the position of this 
peak is assigned the special frequency of the pitch in the data set. This method works well 
30 for low magnification images containing several edge replications in the data set. When few 
repetitions of the periodic structure are present in the data set, the measurement method lacks 
precision. 

[11] While the prior art methods of pitch measurement have proved useful, 
they suffer from several disadvantages. The pitch measurement algorithms that are based on 
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edge detection are more sensitive to noise and hence not very precise in determining the 
pitch. In addition, they require as input to the measurement, the conditions for edge 
detection, including specification of pixel search range intervals for finding each edge. In 
automated pitch measurement, there exists a requirement that the image of the features be 
5 aligned to the search ranges before applying the edge detection algorithm, otherwise the 
measurement will fail. 

[12] Accordingly, there is a need for methods of measuring pitch in the data 
set that can overcome or eliminate such inefficiencies and disadvantages, can be easier to use, 
can result in measurements that are more precise, more robust, more versatile, and do not 
10 require, but can accommodate, search ranges and can be extended to higher data dimensions. 

[13] 

SUMMARY 

4! [14] In accordance with an embodiment of the invention, a method for 

yy measuring pitch in data obtained from metrology and imaging systems is provided. A data 
1=5 set from a metrology or imaging instrument is first obtained. The data set is converted into 
J: digital format if not already in that format. The digitized data set is mapped into a one- 
W dimensional profile data if the digitized data set is not already one-dimensional. The one- 
p dimensional profile data denoted by f(x) is a function of x position values corresponding to 
£\ equally spaced or nearly equally spaced pixels. A criteria function g(T) is constructed as a 
24) one-dimensional data array from the profile data f(x) or any of its derivatives and a 
f7 translation of the profile data f(x) denote by f(x+T) or any of its derivatives. Here, T 
represents the amount of translation, and g(T) is a function of T translation values 
corresponding to equally spaced or nearly equally spaced pixels. A value of translation T is 
then determined either as a whole pixel or with subpixel interpolation such that the magnitude 
25 of g(T) would be either a maximum or a minimum whichever appropriate at said value, 

wherein the determined value is not zero. The determined value is then reported as the pitch 
in the data set. 

[15] In another embodiment, the obtained data set is one or more image(s) 
of portions of a sample, and the one-dimensional profile data f(x) is obtained from the 
30 image(s) through one or more mathematical operations including operations on the obtained 
data set or the one-dimensional profile data to reduce any noise in the data or to shift the 
average of the data or to remove any linear or higher order trend in the data, summing, 
averaging with or without weights, median, and averaging with data culling. 
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[16] In another embodiment, the gray scale image is obtained by one of 
cameras including CCD cameras, optical microscopes, scanning electron microscope 
including top down, tilt, and cross section, scanning ion microscope; transmission electron 
microscope; microscopes in analytical instruments; defect detection and inspection 
microscopes whether optical or scanning electron; microscopes in lithography systems 
including ion beam, x ray, optical, UV, deep UV and extreme UV lithography systems, 
thermal imaging systems, medical imaging devices such as magnetic resonance imaging, 
CAT Scan; ultrasound and other imaging systems such as sonar. 

[17] In another embodiment, the criteria function g(T) is constructed from 
the profile data f(x) and a translation of the profile data f(x) denote by f(x+T) using 
autocorrelation of either the profile data f(x) or any of its derivatives. The computation of 
g(T) is performed with or without normalization of profile data f(x) or its derivatives, with or 
without subtracting a background level from profile data f(x) or its derivatives, with or 
without excluding portions of the profile data f(x) or its derivatives based on their magnitude, 
with or without shifts in the magnitude of the profile data f(x) or its derivatives, with or 
without background uniformity compensation in f(x) or its derivatives, with or without 
allowance for one or more regions in the argument of the profile data or its derivatives that 
can be defined to be excluded from the computations, and by either truncating the summation 
or zero padding or replication to represent the translated values of the profile data or its 
derivatives, or any combination thereof. 

[18] In another embodiment, one or more parameter(s) is reported to 
convey information about the quality of the reported pitch. 

[19] Other advantages and features of the present invention will be apparent 
from the accompanying drawings and the description thereof. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[20] FIG. 1 is a description of prior art method for measuring pitch in the 

data; 

[21] FIG. 2(a) illustrates an example of an image of a portion of a grating 
obtained with a scanning electron microscope; 

[22] FIG. 2(b) depicts the shear and translation operations and is obtained 
from FIG. 2(a) by shearing the top portion of that image and translating it relative to the 
bottom portion by an arbitrary amount; 



# # 

[23] FIG. 2(c) illustrates a special case of shear and translation as described 
in FIG. 2(b) where the translation is equal to the pitch in the image; 

[24] FIG. 3 is a block diagram illustrating steps in a method for calculating 
pitch in the data obtained from metrology or imaging systems in accordance with an 
embodiment of the present invention; 

[25] FIG. 4 is a block diagram illustrating steps in a method in accordance 
with another embodiment of the present invention for calculating pitch in the data obtained 
from metrology or imaging systems; and 

[26] FIG. 5 includes charts comparing the precision of a pitch measurement 
method in accordance with an exemplary embodiment of the present invention and an 
existing pitch measurement method from a commercially available metrology system. 

[27] 

DETAILED DESCRIPTION OF THE INVENTION 
[28] In accordance with an embodiment of the present invention, a method 
for measuring the pitch in any data set obtained from a metrology or imaging system is 
provided where such method does not rely on assigning edge locations, i.e., it does not 
involve application of edge detection algorithms. The method can be extended to data sets 
(images) in any number of dimensions. The data can be obtained from any imaging system 
such as cameras including CCD; optical microscopes; scanning electron microscopes 
including top down, tilt, and cross section; scanning probe microscopes including atomic 
force microscopes and profilometers; scanning ion microscopes; transmission electron 
microscopes; scanning optical microscopes; numerous other microscopes in analytical 
instruments, microscopes in lithography systems; defect detection and inspection 
microscopes whether optical or scanning electron; thermal imaging systems, medical imaging 
devices such as magnetic resonance imaging , CAT Scan; ultrasound and other imaging 
systems such as sonar. 

[29] Broadly, a pitch being defined as the distance between two similar 
objects in the data set, and the data set being represented by f(x) and its translation by an 
amount T represented by f(x+T), the pitch is determined as a specific value of translation 
denoted by P such that the signal f(x) resembles its translated form f(x+P) as closely as 
possible either for a restricted set of values of the argument x or for all x in the data set where 
overlap exists between the original data set and its translation. 

[30] The determination of pitch may include additional conditions such as 
range restriction for P. This definition of pitch is consistent with the period of repetitive 
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structures. Computations without additional analysis or without imposing other conditions 
such as range restrictions can yield values of P that are an integer multiples of the period in 
the data set whenever the original data set contains more than two periods. 

[31] Image shear technique provides a pictorial analogy of what the above 
mathematical construct strives to achieve. Figure 2(a) depicts a top-down image of segments 
of three parallel lines form an imaging instrument. Figure 2(b) depicts the image obtained 
from Fig. 2(a) when the top portion of the image is sheared and translated to the right by 
some arbitrary amount. Figure 2(c) depicts the two images where the translation is equal to 
the pitch in the image. In Fig. 2(c), in the region where the two images overlap, the original 
image and its translation appear almost identical. In image shear techniques, perception by 
human brain is relied upon to measure distances. In accordance with an embodiment of the 
invention, a numerical method for performing a similar task under computer control is 
provided. 

[32] Accordingly, given that a continuous function f(x) represents the 
strength of the signal f as a function of the position x, then the function f(x+T) is the 
translation of that signal by an amount T, a function g(T) referred to as the criteria function is 
constructed from some relationship between f(x) or any of its derivatives and f(x+T) or any 
of its derivatives. The function g(T) is a function of the single variable T denoting the 
translation, and does not depend on any particular value of x denoting the argument of the 
function f( ). For one-dimensional data, an example of the function g(T) can be given as: 

gl (T)= j\f(x)-f(x + T)\dx 

[33] In this example, gl(T) is the integral of the absolute value of the 
difference between the original signal and its translation. Thus, gl(T) is a quantity that 
depends on the amount of translation T, and not any edge position. A particular value of 
translation denoted as P is determined so as to minimize the criteria function, gl in this 
example. That is, gl(P) is the minimum of gl(T) in the range of values considered for the 
translation amount T. 
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[34] Another example for the criteria function g(T) is given below as g2(T) 
which is the integral of the square of the difference between the signal and its translation: 



g 2 (T)= j(f(x)-f(x + T)) 2 dx 



[35] Minimization of the above function results in maximizing the 
autocorrelation function of the original signal f( ). Thus, the position of the peak away from 
10 zero in the autocorrelation function corresponds to the measured pitch value. 

[36] For discrete or digital signal, that is, for a signal f( ) that is sampled at 
regular or nearly regular intervals of its argument x, the above integrals are replaced with 
finite summations. Any criteria functions g(T) will be a discrete function of T. In practice, 
Jpj the pitch P in the data set can be determined to a subpixel level using interpolation. 
fl> [37] FIG. 3 is a block diagram illustrating steps in a method for calculating 

M= pitch in the data obtained from metrology or imaging systems in accordance with an 
?fl embodiment of the present invention. First, in step 40 a data set is obtained from a metrology 
^ or imaging system. Usually, the data sets appear in digital format and are stored in that 
S fashion in computer memory; for example, images obtained from a scanning electron 
f20 microscope or scan lines obtained from a scanning probe microscope. However, if the data 
^ set is not in digital format, it is converted to the discrete form in step 44 where it is 
p represented by an array of data, the array often being one or two dimensional, but possibly 

having more dimensions. If the digital data array is multidimensional, the entire data array or 
a portion thereof is mapped into a one dimensional data profile in step 48. Often straight 
25 averaging is used, as commonly practiced in scanning electron microscopes, to create the data 
profile denoted by f(x) from a portion of the image, but more complicated mathematical 
operations can be used such as averaging with weights. In some circumstances, the data 
array may be aligned prior to averaging to overcome the effects of drift in the data set. The 
data profile f(x) is a discrete function of x position values corresponding to substantially 
30 equally spaced pixels. In step 52, a criteria function g(T) is constructed from the data profile 
f(x) or any of its derivatives and the translation of the data profile f(x) denoted by f(x+T) or 
any of the derivatives of the translated profile. The details of forming the criteria function 
were discussed earlier. Integration or summing over a range of position x values is necessary 
so that the criteria function is not dependent on the value of x as an argument of a function. 
35 In step 56, the location of the peak or trough in the criteria function is determined possibly to 
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a subpixel level using interpolation. The position of the peak or trough, whichever 
appropriate, is reported as the pitch in the data set in step 60. 

[38] FIG. 4 is a block diagram illustrating steps in a method in accordance 
with another embodiment of the present invention for calculating pitch in the data obtained 
from metrology or imaging systems. Step 70 is similar to step 40 of FIG. 3. However, the 
data set is assumed to contain more than a single one-dimensional data profile. Examples of 
such data sets include a sequence of one-dimensional scan lines obtained from a scanning 
probe microscope or a two dimensional image from a scanning electron microscope. Step 74 
is similar to step 44 of FIG. 3 where the data set is converted into digital format if necessary. 
In step 78, the data set is mapped into a number of individual data profiles f(x). In our 
examples, each individual line scan from the scanning probe microscope could be considered 
separately, or the scan lines can be averaged in groups of two or more resulting in a number 
of averaged profiles. A similar construct can be applied to an image from a scanning electron 
microscope where the image is divided in a number of stripes and the scan lines within each 
stripe are averaged to yield a number of one-dimensional data profiles, one for each stripe. 
This is in contrast to the method of FIG. 3 where in step 48, only a single one-dimensional 
data profile was obtained from the data set. For each data profile thusly obtained in step 78, a 
corresponding criteria function is formed in step 82, and the position of either the maximum 
or minimum, whichever appropriate, of each criteria function is found in step 86 in a manner 
similar to step 56 of FIG. 3. Finally, in step 90, the pitch in the data set is reported as either 
the average of all such positions, or another suitable statistical function of such positions, for 
example, median. 

[39] The methods of present invention cover maximization or minimization 
of any reasonable criteria function g(T) to determine pitch. This approach avoids individual 
edge position assignment and instead, uses a range in x of the signal f( ) to determine the 
pitch. The present invention includes but is not limited to maximum autocorrelation of the 
signal. Other examples are autocorrelation of the derivative or functions similar to gl as 
mentioned before, or any of these functions with weights included in the integrand. For any 
given criteria function g(T), any number of methods can be used to find the maximum or the 
minimum. As an example, when the criteria function is autocorrelation, Fourier transform 
methods together with application of ParsevaPs theorem can be used to find the 
autocorrelation peak position in space from the position of the peak in spatial frequency 
domain. 
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[40] Maximum cross correlation technique is commonly used in several 
disciplines such as time delay measurement, overlay measurement, and pattern recognition, 
that is, finding objects in a target image. Based on the known properties of cross correlation 
and depending on the nature of the detected signal, the autocorrelation technique can be the 
method that provides the lower bound on the precision of determining the pitch in the data 
set. That is, the autocorrelation algorithm can produce pitch measurements with best 
precision for a given set of experimental conditions. An example of improved pitch 
measurement precision is given in Figure 5 where the results are compared between a 
measurement algorithm of a commercially available scanning electron microscope and an 
exemplary embodiment of the present invention. The experimental conditions were typical 
of metrology applications in semiconductor manufacturing industry. In arriving at the results, 
an identical data set was used for both measurement methods. The metrology tool reports the 
results of measurement and also stores in memory the profile data used to arrive at the 
measurement. The stored data profiles for all the measurements were retrieved and analyzed 
off-line with the method of present invention. Eight different replications of the same feature 
on a silicon wafer were measured. The wafer was loaded into the system, measured, and 
unloaded from the system twice using an automated recipe-driven operation. In each cycle, 
each of the eight features was measured 8 times consecutively, where no feature was moved 
in between the eight consecutive measurements. For the commercially available 
measurement method, the pitch measurement precision was slightly more than 1% of nominal 
size (3 standard deviations). This value is consistent with the precision of pitch measurement 
for similar metrology tools as reported in the literature in which the 3 standard deviation of 
measurement ranges from 0.6% to 2% of measured value. The range of pitch measurement 
precision values can be found in an article entitled "Evaluation of the long-term stability of 
critical-dimension measurement scanning electron microscope using a calibration standard" 
by F. Mizuno et al, Journal of Vacuum Science and Technology B, Vol. 15, No. 6, pp. 2177- 
2180, 1997. For the present experiment, the results for the exemplary embodiment of the 
invention were improved compared to the commercially available measurement method by a 
factor of 3. In this exemplary embodiment, the criteria function was obtained by normalized 
autocorrelation of profile data with background subtraction to result in zero average value of 
the profile data. 4 

[41] Though improvement in precision is one attribute of the present 
invention, it is not limited in its application to that purpose. Another issue is the ease of 
measurement and the requirements for the measurement setup. Prior art measurement 
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methods require defining intervals in the argument x (restricted areas of the image or data set) 
within which the edge of interest must be found. The methods of present invention, while 
allowing for this option, do not make it a particular requirement. Thus the methods can be 
modified to construct the criteria function from the ranges of interest, and hence be a drop in 
replacement for the operation of existing pitch measurement algorithms. 

[42] Furthermore, the methods of present invention do not rely on the 
sharpness of the image of the edge. When prior art techniques are used for magnification 
calibration, samples with sharp physical edges are preferred and often required by system 
manufacturers. Thus, magnification standards are currently required to possess nearly 
vertical sidewalls so that their cross sections resemble trains of rectangular pulses. In 
accordance with the present invention, the need for samples with nearly vertical edges for 
magnification calibration is alleviated. The methods of the present invention would work 
with periodic structures of any cross section shape such as triangular pulses. In addition, in a 
light microscope image at high magnifications approaching the resolution limit of light optics 
where the edge profile is spread due to diffraction and other limitations of a microscope, the 
pitch measurement can still be easily performed with the methods of the present invention. 

[43] The methods of the present invention can further be modified to 
include only certain portions of the profile data, based on the magnitude of f(x), or 
incorporate shifts in the signal intensity to add or remove bias from the signal. The criteria 
function may be computed with or without normalization of the profile data, with or without 
background subtraction, with or without background signal uniformity compensation, with or 
without allowance for regions in x of the one-dimensional profile data f(x) that can be 
defined to be excluded from the computations, and by either truncating the summation or 
zero padding or replication to represent the translated values of the one-dimensional profile 
data. 

[44] The method of the present invention can easily be extended to measure 
pitch in the data set in two or more dimensions. In addition to one-dimensional gratings, 
samples containing two-dimensional arrays of equally sized squares or equally sized disks 
can be used in metrology and imaging systems. Arrays of squares can be used for 
establishing the position of the scanned electron beam relative to the stage in electron beam 
lithography systems. Arrays of disks are used in calibrating electronic circuit test systems 
that employ ball grid arrays. The methods of present invention, as discussed above, can be 
applied to each dimension of the data set to obtain the pitch in the data set in that dimension. 
Alternatively, the criteria function can be constructed in two or more dimensions and the 
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search for the position of the maximum or the minimum of the criteria function can be 
conducted in two or more dimensions. 

[45] While the present invention has been particularly shown and described 
with reference to the preferred embodiments, it will be understood that the various changes 
and modifications will occur to those skilled in the art without departing from the scope and 
true spirit of the invention. The scope of the invention is therefore to be determined solely by 
the appended claims. 
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